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Orthogonal frequency division multiplexing (OFDM) is a multicarrier 
modulation (MCM) technique that divides the wide bandwidth into parallel 
narrow bands, each of which is modulated by orthogonal subcarriers. 
Currently, OFDM is a high-spectral efficiency modulation technique that is 
used in a variety of wired and wireless applications. The transmitted signal 
in a wireless communication channel spreads from transmitter to receiver 
through multiple reflective paths. This triggers multipath fading, which 
causes variations in the received signal's amplitude and phase. Slow/fast and 
frequency-selective/frequency-nonselective are the main types of multipath 
fading channels. Therefore, in this paper, we proposed new models for 
modeling multipath fading channels, such as the exponential fading channel 
and the Gamma fading channel. In addition, new bit-error-rate (BER) 
derivations have been derived. The performance of the OFDM system over 
proposed channel models has been evaluated using Monte-Carlo simulation 
and compared to the Rayleigh fading channel model. The obtained results 


via simulations show that the exponential fading channel at a rate parameter 
(A=0.5) outperforms the Rayleigh fading channel by 6 dB for all values of 
Eb/No, while the Gamma fading channel at (a=2) outperforms the Rayleigh 
fading channel by 3 dB for all values of Eb/No. 
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1. INTRODUCTION 

The growing interest in wireless broadband communications and interactive media has fueled 
intense study on data rate transmission, bandwidth efficiency, and bit-error-rate (BER) performance across 
multipath fading channels [1]. Since the information to be sent is modulated into a single carrier in a 
conventional communication system, a single fade or interferer may cause the whole connection to fail. To 
achieve high bit rates, the symbols must be sent rapidly, using the full bandwidth. When a channel is 
frequency selective, its impulse response may span several symbol periods, resulting in what is referred to as 
intersymbol interference (ISI). This is a significant issue in wideband transmission across multipath fading 
channels, since it makes it more difficult to identify the sent signal, thus degrading the system's performance 
significantly [2]. Orthogonal frequency division multiplexing (OFDM) has been explored for a variety of 
wireless (applications, standards) and wired because of its popular features of high bandwidth efficiency, 
easy channel equalization, and immunity to multipath fading channels [3], [4]. This modulation method may 
be utilized effectively and consistently for high-speed digital data transmission across multipath fading 
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channels, since it splits the available bandwidth into multiple orthogonal data streams [5]. As a result, fading 
will impact a very small proportion of the subcarriers. Each subchannel will occupy a small frequency 
interval in the frequency domain, where the channel frequency response will be almost constant and each 
symbol will experience a nearly flat fading channel. Therefore, it is seen as one of the most promising 
solutions to the ISI problem [2]. Intercarrier interference (ICT) is eliminated when orthogonal subcarriers are 
used, and the signals carried on the various subcarriers do not interfere [6]. Hence, a simple frequency 
domain equalization method can be used [7]. IEEE 802.11 a, g, and n, IEEE 802.16a, and terrestrial digital 
video broadcasting (DVB-T) have all adapted it [8] and ODFM is used in digital subscriber line (DSL) and 
asymmetrical digital subscriber line (ADSL) for wired applications [9]. In the last few years, there have been 
a lot of studies done on the success of wireless channel models. In [10], the study of the performance 
evaluation of the OFDM-based 256 and 1024-QAM under multipath fading propagation circumstances 
defines the acceptable signal-to-noise ratio (SNR). In [11], focuses on the performance evaluation and 
investigation of OFDM systems operating across multipath fading channels like additive white Gaussian 
noise (AWGN), Rician, and Rayleigh. In [12] OFDM is studied and simulated in the presence of Rayleigh 
fading channels using channel estimation and equalization. Bit error probability is used to quantify 
performance. In [13], it analyzes the bit error rate (BER) performance of discrete cosine transform (DCT)- 
based OFDM and compares it to discrete fourier transform (DFT)-based OFDM to determine which 
modulation scheme performs better. In [14] The cyclic prefix's effect on the OFDM system's performance 
was investigated. The BER for different lengths of the cyclic prefix was compared using an OFDM 
transceiver equipped with a 16-QAM modulation method, highlighting the effect of channel noise. 

In this paper we suggested new models to find new simulations and effects as well as to improve the 
BER performance of the OFDM system. In addition, two fading channel scenarios: exponential fading channel 
and Gamma fading channel have been adopted in this paper. The major contributions to this paper are related 
to improving the signal detection and enhancing the BER performance analysis of 4-QAM/OFDM over 
Gamma fading channels and exponential fading channels in AWGN based on the following: i) we derive the 
real part, the imaginary part of the channel impulse response in the time domain, in addition to the effective 
probability distribution and the effective SNR distribution of the Exponential and Gamma fading channels in 
the frequency domain at the OFDM demodulator output by using the central limit theorem (CLT), ii) we 
derive the variance in each dimension, i.e., the real part and the imaginary part of the exponential and 
Gamma fading channels, iii) we derive the BER based on the effective derived SNR distribution in the 
Gaussian noise only and in the presence of the exponential and Gamma fading channels in AWGN, and iv) 
we detect the received signal by using an optimum detector based on the Euclidean distance depending on the 
derived effective probability distribution at the OFDM demodulator output. 

The rest of this paper is organized as follows: section 2 system model, while section 3 is the 
derivation of the variance in the time domain for exponential distribution and Gamma distribution. Section 4 
BER performance analysis of M-ARY digital modulation over AWGN, exponential and Gamma fading 
channels, section 5 describes the simulation results and section 6 concludes this paper. 


2. RESEARCH METHOD 
2.1. System model 

In this paper, the binary data bits dẹ are divided into two-bit classes before being converted into the 
2? symbols of a 4-QAM constellation to obtain the modulated symbols X=M(c) according to constellation 
mapping c shown in Figure 1 [15]. The complex base-band OFDM signal x„ can then be obtained in the time 
domain using an N-points inverse fast fourier transform (IFFT), as seen in [16]. 


pty eN, 0<n<N-1 (1) 
n VN i=0 ^l , = = 

where X; the complex modulated symbol and N is is the number of orthogonal subcarriers. 

To minimize ISI between consecutive OFDM symbols in multipath fading channels, a time-domain 
CP of duration Ncp samples is added at the beginning of each OFDM symbol to produce the transmitted 
symbol £ = [xy — Ncp, Xy — Ncp + 1, ..., Xy — 1, Xo X1 Xy — 1], which is programmed to exceed the 
cumulative delay spread of the fading channel. Then, ¥ is transmitted through the fading channel [17]. This signal 
is combined with additive noise in the wireless channel, resulting in the obtained signal as shown in (2) [18]: 


y = h@k+w (2) 


where ® represents the convolution operation, h represents the complex impulse response coefficients in the 
time domain of the multipath fading channel and w is the AWGN. 
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Figure 1. The four signal points of 4-QAM modulation in the constellation diagram 


The multipath fading channel of wireless communication was modeled using exponential 
distribution or Gamma distribution in [19], [20] due to the fact that a large number of distributions used to 
model the multipath fading and shadowing fading channels for survival analysis are special instances of the 
generalized gamma distribution, such as the exponential and the gamma distributions [19], [21], as explained 
in more detail in this section. If the amplitude of the complex random variables h of the channel impulse 
response follows the exponential distribution in the time domain, the distribution can be expressed in [20] as 
in (3): 


fah) = Ae" 0 <h<o, A>0 (3) 


where 4 is the rate parameter of the fading channel. 

The main concern with the fading can be stated into real and imaginary parts expressed as h, and h;, 
respectively. Assuming the magnitude of the fading channel h in the time domain, |h| = y h2 + h?, exhibits 
an exponential distribution where h = h, + jh;. Two random parameters are expanded by 0, which is a 
randomly distributed random variable between the limits of the bounds z and -z, creating two further random 
variables as in (4) and (5). 


h, = |h| cos(@) (4) 
h; = |h| sin(@) (5) 


As a result, the PDF of the real part of the multipath fading channel can be expressed as (6) [22]. 


Ir cos(0) (6) 


dh 
Then, f (h,)|@ a conditional PDF of h, subjected to 0, can be written as (7) [22]. 


—f™®_ _ 4,-an__1 
fh, )|8 = dh,/dh de dh,/dh 
-àh 


ecos(®) (7) 


= jeos(6)| 


Moreover, f (h;)|@ a conditional PDF of h; subjected to 8, can be written as (8): 


Ee a (ee ee | 
f (hi) = dhi/dh =4e dhi/dh 
—Ah 


pare (8) 


= jsin(6)| 
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Then the complex received signal, y, processed using the CP removal and fast fourier transform (FFT) 
operation [18]. As a result, in the frequency domain, the complex signal can be represented as (9): 


Y=HX+W (9) 
Y, + jY; = (H, + jH;)(X; + JX) + (W, + jW;) 


where Y, H, X, and W represent the complex received signal in the frequency domain, the transfer function 
of the complex fading channel, complex modulated signals in the frequency domain, and AWGN, 
respectively. 

Hence, the distribution of the real part, f(H,) and the imaginary part, f(H;) would approach a 
Gaussian distribution after the FFT operation with a mean equal to zero, u=0, and variance Oi, and Oh 
computed in the next section. According to the central limit theorem (CLT) [23] as (10). 


Hy? 
fH) = -e r (10) 
200k, 
H;? 
f(H) =——e i (11) 
2h, 


As a result, the magnitude of the fading channel in the frequency domain can be computed as |H| = 
y H” + H,? would follow the Rayleigh distribution with variance ø computed in the next section as (12): 


|H|2 


fH) = He 2h (12) 
on 


with uniformly phase distributed in the range [-7, 2). Moreover, if the amplitude of the complex random 


variables h, |h| = ~h? + h? of the channel impulse response follows the Gamma distribution in the time 
domain, the PDF can be expressed in [20] as (13): 


1 


h 
a-la p 
op e F0<h<w,a>0,8>0 (13) 


falh) = 


where a is the shape parameter and ß is the scale parameter. Following the same mentioned discussions, the 
distribution of f (h,)|@ and f (h;)|@ can be derived as (14) and (15). 


h 
- O a-1 o`p 
f(h,)18 = dhy/dh  T(a@)p% Beene dh,-/dh 
1 n yoI hk 
SEER ft Bcos(@) 
T'(a)B™|cos()| (5) 4 ag 
and 
f(hi)jo = LE = — ho E E S 
dhi/dh  T(a)R% dhi/dh 
1 h \o1 FETO ao 
z= ——— | — sin 
T(a)B%|sin(@)| (5) e (15) 


Thanks to the central limit theorem after FFT operation, which makes the distribution of the real and 
imaginary components of the Gamma fading channel approaches Gaussian distribution as expressed in (10) 
and (11) in the next section with u=0 and variance of, and On, computed in the next section. Therefore, the 
magnitude of the fading channel with Gamma distribution will follow Rayleigh distribution as given in (12) 
and the phase will follow uniform distribution throughout the frequency domain from the range [-7, 7). 
Finally, the received signal in (9) can be detected using maximum likelihood (ML) detection. The ML can be 
implemented by finding the minimum Euclidean distance between all the possible transmitted symbols over 
the channel H and the received signal Y as [24]: 
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where C=C1, C2, C3, C4 represents a collection of all possible symbols X in a 4-QAM constellation as 
presented in Figure 1. 

The bit error rate (BER) is determined by comparing the binary data bits generated by the source at 
the transmitter, dk, with the eventually obtained bits at the receiver, dp, and then dividing the result by the 
number length of dk [25]. Figure 2 illustrates the configuration of an OFDM transceiver system (transmitter 
and receiver) [26], [27]. 
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Figure 2. Baseband OFDM system model over exponential/Gamma fading channel 


2.2. Derivation of the variance in the time domain for exponential distribution and gamma distribution 
2.2.1. Variance derivation of exponential distribution 

Suppose the magnitude of the complex random variable h in the time domain follows the 
exponential distribution, that is, it has pdf given by (3). Hence, the variance can be computed as in (17) [20]: 


on = E(\h|?) — |E(h)|? (17) 


where E(|h|?) is the second moment and |E (A)|? is the square of the first moment (expected value). 

The variance of an exponential distribution is computed in [20] as of = = However, the mean and 
variance of a complex random variable, h = h, + jh;, can be expressed as given in [24] as E(h) = E(h,) + 
jE(h;) and of = Oi, + Oh, So, the variance of the real component Oi,» is equal to the variance of the 


imaginary component Oi, and the variance in each dimension is: 


z 1 
Oi, = of, == (18) 
2.2.2. Variance derivation of gamma distribution 
If the magnitude of the complex random variable h follows the Gamma distribution as given in (13), 
the variance can be expressed as given in [28] as of = af, hence, the variance in each dimension can be 
expressed as mentioned above as (19). 
2 2 
On. = Of, =t- £ (19) 
2.3. BER performance analysis of QAM over AWGN, exponential and gamma fading channels 
The influence of fading is measured using 4-QAM modulation. The bit error rate (BER), is a more 
accurate performance parameter for measuring modulation schemes. In a slow flat fading path, the following 
integral may be used to measure the BER performance of any digital modulation scheme [29]. 


Pp= h Pp awan (y) Par (y) dy (20) 


The error probability of a specific modulation scheme in an AWGN channel at a specified signal-to- 
noise ratio is: 
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y= (21) 


No 


where h is the random variable of channel gain, in a non-fading AWGN channel, sb is the ratio of the bit 


energy to the noise power density. The random variable h? denotes the fading channel's instantaneous power 
and Par(y) is the probability density function of y due to the fading channel [29]. 


2.3.1. Derivation of BER of 4-QAM digital modulation in AWGN channel 


The cumulative distribution function (CDF) of a continuous random variable X can be expressed as 
the integral of its probability density function py follows [20]: 


Fy(x) = f°, px (edt (22) 


The probability of error receiving the constellation symbol C2=1+1j over AWGN channels as shown 
in Figure 1 can be computed as (23) [30]. 


P;(C2) = F¥?(0) + Fy?(0) — Fy?(0)F¥?(0) (23) 


Because of the symmetry, the marginal probabilities for the real and imaginary components can be calculated 
using the Gaussian PDF as (24). 


Fy? (0) = Fy?(0) (24) 


The equation of Gaussian PDF is: 


p(x) = age, with u = 0 and o? = xe (25) 
So 

F£? (0) = tee — e 2% dx= tert =| =erf( 00) 

F£? (0) = z- tert( fE), where y = ee (26) 


As shown in (26) can be rewritten as (27): 


FẸ) = Lerfe( B 27) 


Now we substitute (27) into (23) we get: 


EEEE 
=erfc (J -terfe?( z) 


= erfc (H fory>0 (28) 


Due to the constellation symmetry in Figure 1, the error probability of C2 is equal to the error probability of 
CO, C1 and C3. The total symbol error probability can be expressed as (29). 


Ps(C) = P (CO) Pe (CO) + P (C1) Pe (C1) + P (C2) Pe (C2) + P (C3) Pe (C3) (29) 


We assume that four symbols are transmitted with equal probability, so P(CO) = P(C1) = P(C2) = P(C3) = 
: , hence, the net symbol error probability can be expressed as (30). 
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Ps(C) = =(Pe(CO) + Pe(C1) + Pe(C2) + Pe(C3)) = + (4erfe (E) ) 


=erfc (E) (30) 


Therefore, the BER of the 4-QAM constellation can be computed as (31): 


— Ps(C) _ Ps(C) 
Pb(C) = log2M pe a 
Pb(C) = serfc( |22- ) (31) 


Leads to the 


1 E 
BERggam,awen =3erfe ( B (32) 


2.3.2. Derivation of BER of 4-QAM digital modulation in exponential fading channel 

For the exponential fading channel as expressed in (3) in the time domain, the H follows the 
Rayleigh distribution in the frequency domain after FFT operation based on CLT as computed in (12), and 
|H]? has a chi-square distribution of two degrees of freedom as expressed in (33) [29]: 


Pary) =5exp( -7 ) (33) 


where 7 = = E(|HI?] is the signal-to-noise ratio on average for E[|H|*] = 1/42, substituting in (33), gets 


“F 


zoel GIRS Epia 
Par) = GDE ° aie ~~ 


Substituting (32) and e (34) in (20) we get: 


on 7 SO BSS 
Pp= de serfc(y) * GATS e G)Ep/No dy 
y(1/A?)(Ep/No)+1) 
(1/2) Ep /No TE 


eee! Gee 
P _1 Gz S)Ep/No | 
2 JA) (Ep /No) +1 serfo) e 


0 


The equation after substituting the integral values becomes: 


0 
6) [Ep/No erf(o0) 1 * er fc(oo) $ a om (1/A)VEb/No erf(0) 1 er fe(0)e T27?" 
2/(1/42)(Ep/No)+1 2 
2 JE Ep/No)+1 
where erf (œ) = 1,erfc(o) =0,erf(0) = 0,erfc(0) =1 
_G/AVE5/No__ 1 (35) 


© 24 G/A2)Ep/No)+1 2 
As shown in (35) can be rewritten as (36). 


-ifj EM 
m= ( ee) (36) 


2 


2.3.3. Derivation of BER of 4-QAM digital modulation in gamma fading channel 
Using the same steps that are used to find the Paf (y) of the exponential fading channel, the Pas (y) 
for the Gamma fading channel is: 
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>Y: 


Pap) = earns (37) 


1 
aB?E,/No 


Sub (32) and (37) in (20) we get: 


-Y 
œ 1 1 ZBZFE, /N 
Prha TN * am y 


{V(@B? (Ep/No)+1) 
VaB2*Ep/No erf| ——————$——— 
b/No í aß2»*Ep/No j 


-Y 
=| —— m Serfe(Vy) e %B7+(Ep/No) H 


2 @B2(Ep/No)+1 


The equation after substituting the integral values becomes: 


= | yaß?*Ep/No erf()  erfc(Vo) 7 aB2*(Ep/No) | yaß?*Ep/No erf(0) 
2/aß?(Ep/No)+1 2 2 @B2(Ep/No)+1 


0 
~erfic(v0)e%**p/No) | 
V @B2*Ep/No 1 (38) 


2 (@B2*Ep/No)+1 2 


As shown in (38) can be rewritten as (39). 


_1f, _ VaB?*Ep/No 
Po =] (1 - a] 69 


3. SIMULATION RESULTS AND DISCUSSION 

The BER performances are investigated using the MATLAB simulation program for 
4-QAM/OFDM over two different channel distributions such as the multipath exponential fading channel and 
the multipath Gamma fading channel in AWGN, and the results are compared with the 4-QAM/OFDM over 
the Rayleigh fading channel. The following simulation parameters were set, the number of sub-carriers was 
set as N=64 and CP=16, the constellation size was set as 4-QAM and the system bandwidth was set at 
20 MHz for the ITU pedestrian standard channel for multipath exponential fading channel and multipath 
Gamma fading channel are used, as shown in Table 1 [31]. 


Table 1. ITU channel profile 
ITU Pedestrian Ch.103 
Pathn Power (dB) Delay (us) 


1 0 0 

2 -0.9 0.2 
3 -4.9 0.8 
4 -8 1.2 
5 -7.8 2.3 
6 -23.9 3.7 


3.1. Performance of 4-QAM/OFDM system over exponential fading channel 

We are interested in this section to show the matching of the derived distributions with simulated 
results in the time domain based on (3) and the frequency domain based on e (12), for oå computed in (18). 
In Figure 3, we see the exponential fading histograms, where the parameters are set to A=0.5, A=1 and A=1.5, 
respectively. It should be noted that (3) and the derived closed-form PDF in (12) have similar analytical 
distributions to those of the simulation results in the probability domain. 

The impact of changing à parameter on the BER performance of the 4-QAM/OFDM system over 
exponential fading channels and AWGN compared to the Rayleigh fading channel is shown in Figure 4. 
Simulations were run over different scenarios of the channel effect by changing the scale parameter À to 0.5, 
1, and 1.5 for severe impact, normal impact, and low impact, respectively. For instance, in the case of the 


Int J Elec & Comp Eng, Vol. 12, No. 5, October 2022: 5123-5135 


Int J Elec & Comp Eng ISSN: 2088-8708 O 5131 


normal impact, the exponential fading channel has the same effect as the Rayleigh fading channel, and in this 
case, the parameter of À is chosen as follows: 

The PDF of the signal to noise ratio distribution at the OFDM demodulator output can be derived 
based on CLT as given in (34), while in the Rayleigh fading channel it can be expressed as in [32]. 


-Y 
= 1 (202) Ep /N 
Pas) ease biNo (40) 


Hence, both distributions have the same effect when 
= 207 
By take the square root for both sides: 


ul 
{7 V20 
a 


The standard deviation o = 
v2 


in both directions of the Rayleigh distribution. Hence, 


i = V2 T soda=1 

Based on the first value (A =0.5), we determined that the BER performance of the 4-QAM/OFDM 
method over an exponential fading channel is superior to that of the system over a Rayleigh fading channel. 
At the second value (A =1), the BER performance of the exponential fading channel is close to or similar to 
that of the Rayleigh fading channel. But on the third value (A =1.5) the performance will be worse. Table 2 
shows the BER performance of the proposed system compared to the conventional system at Eb/No=20 dB. 
The Table 2 shows that the BER improved by decreasing the rate parameter à, and the performance is 
matched to the conventional system over Rayleigh fading channel when A=1. 


Exponential distribution in the time domain 
T T T T T T T 


Exponential distribution in the frequency domain 
T T T T T T T 


Figure 3. Simulated Vs histogram plots in the time domain and the frequency domain for 
exponential fading channel 
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Table 2. BER values with different scenarios of à parameter 


à parameter BER (Exponential) BER (Rayleigh) 
0.5 0.6238 x10° 0.2481x107 
1 0.2481x107 0.248110 
1.5 0.5532x107 0.2481x107 


3.2. Performance of 4-QAM/OFDM system over gamma fading channel 


A comparison 


of the histogram plots of the derived distributions with simulated results for Gamma 


distribution as presented in the time domain in (13) and in the frequency domain will follow the Rayleigh 
distribution as given in (12), for of computed in (19). In Figure 5, we see the Gamma fading histograms, 
where the parameters are set to o=0.5, a=1 and a=2 for B=1. It should be noted that (13) and the derived 


closed-form PDF in ( 
probability domain. 


12) have similar analytical distributions to those of the simulation results in the 


3 Gamma distribution in the time domain 
T T T T T T 
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Figure 5. Simulated Vs histogram plots in the time domain and the frequency domain for 
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The impact of changing a and P parameters on the BER performance of 4-QAM/OFDM systems 
over Gamma fading channel and AWGN compared to Rayleigh fading channel is shown in Figure 6. 
Simulations were run for various values of the channel effect by changing the a parameter (0.5, 1, 2) for 
severe impact, normal impact and low impact, respectively, and for the B parameter equal to 1, by using the 
same method to find the normal effect in the exponential fading channel, we found the normal effect in the 
Gamma fading channel at o=1. 
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Figure 6. Performance of 4-QAM/OFDM system over gamma fading channel 


In this figure, the BER performance of the 4-QAM/OFDM system over a Gamma fading channel is 
lower than that of the system over a Rayleigh fading channel based on the first value (a=0.5). Moreover, the 
BER performance of the Gamma fading channel is equal to or equivalent to that of the Rayleigh fading 
channel at the second value (a=1). Furthermore, the result will be better on the third value (a=2). Table 3 
shows the BER performance of the proposed system compared to the conventional system at Eb/No=20 dB. 
The Table 3 shows that the BER improved by increasing the shape parameter a, and the performance is 
matched to the conventional system over Rayleigh fading channel when a=1. 


Table 3. BER values with different scenarios of Shape parameter (a) at B=1 
a parameter BER (Gamma) BER (Rayleigh) 


0.5 0.4899 x10? 0.2481x102 
1 0.2481x102 0.2481x102 
2 0.1245x107 0.248110? 


4. CONCLUSION 

The BER of the 4-QAM/OFDM system in AWGN over different channel scenarios, such as 
exponential and Gamma fading channels, in addition to the BER performance, has been investigated for both 
channels in this paper. The simulation results obtained via the MATLAB program have confirmed the 
feasibility of the exact derived PDFs at the OFDM demodulator output in the frequency domain based on 
CLT, leading to derived exact new theoretical BER equations matching to the Monte-Carlo simulation for 
both proposed channel distributions. The 4-QAM OFDM system in AWGN over exponential fading channel 
with different à parameters (0.5, 1, and 1.5) has been simulated. It has been shown from the results that the 
exponential fading channel outperforms the Rayleigh fading channel when (A<1) for all values of Eb/No. 
Additionally, the 4-QAM OFDM system in AWGN over Gamma fading channel with different values for a 
(0.5, 1, and 2) and for B=1 has been simulated. It has been shown from the results that the Gamma fading 
channel outperforms the Rayleigh fading channel when (a >1) for all values of Eb/No. 
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